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Abstract 

Application of the path-integral approach to continuous measurements leads 
to effective Lagrangians or Hamiltonians in which the effect of the measure- 
ment is taken into account through an imaginary term. We apply these con- 
siderations to nonlinear oscillators with use of numerical computations to 
evaluate quantum limitations for monitoring position in such a class of sys- 
tems. 
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Several numerical experiments have shown that quantum-classical correspondence in 
chaotic dynamics seems to be anomalous [1]. In particular, it has been suggested that 
nonlinear systems showing chaotic behaviour in a classical regime may exhibit suppression 
of their stochasticity when quantization is taken into account [2]. However this suppression 
should appear on a time scale so short that the classical chaotic behaviour could not be 
observed in real systems. To recover the Liouville dynamics use of quantum measurements 
process has been suggested [3]. In [4] an example of a chaotic system in a quantum regime 
under instantaneous measurements has been investigated. 

The quantum measurement processes for continuous monitoring of some observables of 
a dynamical system may be taken into account through path-integral formalism [5-7] and 
offers the possibility to investigate in a general way suppression of quantum behaviour. It 
is therefore interesting to apply this method to nonlinear systems, specifically those pos- 
sessing chaotic properties in the classical regime, having as main goal the study of quantum 
suppression of chaotic behaviour in such a class of systems. In this letter we deal with the 
development of a technique based upon numerical integration of an effective Schrodinger 
equation and applicable to continuous quantum measurements in nonlinear systems. As an 
example a nonlinear oscillator will be investigated. 

The path-integral approach to continuous measurements is essentially based upon re- 
striction of Feynman path integrals [5]. The integration is restricted to the set of paths 
compatible with the informations, namely the given output of the continuous measurement. 
If the system moves between the points of the space-time (x',0) and (x",t) the result of 
integration K\ a ](x", r; x', 0) depends on the measurement output [a] = {a(t) | < t < r}. 
This can be interpreted in two different ways. Firstly, it represents the probability amplitude 
for the measurement output a(t) in the time interval r given the positions of the system 
x' and x" before and after the measurement. Secondly, it can be seen as a propagator for 
the system subject to the continuous measurement given the output of the measurement. 
The restriction of the path-integral on some set of paths compatible with the measurement 
output can be effectively done with the help of a weight functional W[ a ] [x] depending on the 
measurement output a(t) and decaying outside the set of paths compatible with a(t) 

K [a] (x",T;x',0) = J d[x]exp^J L(x,x,t)dt^w [a] [x] (1) 

For example, if the coordinate x is monitored with error Aa and the result of the measure- 
ment is a(t), the weight functional W[ a ][#] selects those paths x(t) in the corridor centered 
around a(t) and having a width Aa (see Fig. 1). A simple choice is a Gaussian weight 
functional 

w [a] [x] = exp{-^^ J\x{t) - a(t)) 2 dtj (2) 

The restricted path integral is then 

K[ a ](x",T;x', 0) = J d[x]exp^J^ L(x,x,t)dt - (x - a) 2 dtj (3) 

It is important to observe that the resulting path integral may be considered as describing 
a free (i.e. not measured) system but with effective Lagrangian having an imaginary term 
due to the measurement [8] 
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L eff (x, x, t) = L(x, x, t) + —<^( x ~ a^)) 2 ( 4 ) 

This contribution produces a decrease of the density of the system in the configuration space 
far from x(t) = a(t). The decrease is linked with the restriction of the alternatives due to 
the measurement performed [9]. Finally the effective Lagrangian is time-dependent even if 
the original Lagrangian is not. 

In order to estimate a probability distribution for the measurement output of a physical 
system we consider the convolution 

I[a] = ((h\K[a]\<f>i) = J J (f>* 2 (x")K [a] (x",T;x> ' ^^(x^dx'dx" (5) 

According to the first interpretation of K\ a ] the quantity Jui is a probability amplitude for 
the measurement to give the output a(t) under the condition that the system has been in 
the state <f>i before the measurement and in the state 2 after a time r. The probability 
distribution for the measurement output is then 
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Note that P[ ] depends upon the instrumental uncertainty Aa. 
The amplitude J[ ] can be also written as a scalar product 

I [a] = toM«](r)> (7) 

where 

i> [a] (x",r) = jK [a] (x",T;x',0)Mx'W (8) 

According to the second interpretation of K[ a ] the wave function ip^ (x, t) represents the 
evolution at time t of the state 4>\{x) under the action of a continuous measurement with 
output a(t). This also means that ip[ a ](x,t) can be found as the solution of the time- 
dependent Schrodinger equation 

in d -^^ = H eff ^ [a] {x ) t) (9) 

with an effective Hamiltonian H e ff corresponding to the Lagrangian L e ff in (4) and with 
the choice ^\ a ](x, 0) = <f>i(x). 

The time dependence of the effective Hamiltonian H e ff makes analytical calculations dif- 
ficult. In this case one may more simply evaluate Pr a i through the Feynman propagator K\ a y 
However, due to the quadratic nature of the measurement contribution to L e ff, analytical 
calculations are essentially restricted to the case in which L(x,x,t) is a linear oscillator. In 
general a numerical approach must be followed and in this case the Schrodinger formalism 
is more suitable. The partial differential equation (9) is reduced to a simple finite difference 
recursive equation by choosing a proper lattice to simulate the continuous space-time [10]. 
No particular problems arise from the time dependence and the non-Hermitian nature of 
the differential operator H e ff [11]. 



3 



In order to understand the behaviour of a nonlinear system let us first consider a linear 
oscillator 



m . 2 moo 



In this case analytical results have been obtained allowing a test of the numerical technique. 
The effective Lagrangian corresponds to a forced linear oscillator 



m . 2 mu 2 2 2ih ih 



with renormalized complex frequency 

«Z = - 2 - (12) 

mrAa z 

For any choice of the measurement output a(t) the propagating kernel Kui can be easily cal- 
culated [12]. Let us consider what happens when the quantum system is in the ground state 
of the unmeasured oscillator before and after the period r of the continuous measurement 

AM = *,(,)=(_) exp(-— *») (13) 

Due to the shape of 0i and 02 it is natural to choose the null boundary conditions a(0) = 
a(r) = for the measurement output a(t). Any such function a(t) can be written as a 
Fourier sine series. We consider only measurement outputs of the form 

71 

a(t) = e sin Q n t ; fl n = n— (14) 

T 

where n is an integer number. For a fixed Q n the probability distribution P[ ] is reduced to 
a function of the amplitude e. P(e) is a Gaussian function of width Aa e // 

1 / e 2 \ 

P ( e ) = ^^ exp ( — - — ^ — ) (15) 

where 

AaZh = 2 3?e 
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muuT 2 Aa\ttl -uj 2 ) 2 



, .^r/ , x ("I)" 

1 — « — cot(a; r r) + 



V sin(ccvr) 



(16) 



The meaning of Aa e // is linked to the role of the measurement device during the evolution 
of the system under measurement. When the instrument error Aa is large in comparison to 
the characteristic quantum scale of the system under measurement, ^Jh/mcu in this case, the 
corridor shown in Fig. 1 contains the classical trajectory (with null boundary conditions). 
Thus the classical limit is 
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Aa c = lim Aa e // = Aa (17) 

Aa-^oo 

On the other hand when Aa becomes small quantum noise arises. Also corridors far from 
the classical trajectory are probable. The quantum limit is written as 



-1/2 



(18) 



i.e. the effective error Aa e jj diverges as Aa~ 1//2 . In an intermediate situation Aa e // interpo- 
lates between these two limits always maintaining values larger than the instrumental error. 
In Fig. 2 the behaviour of Aa e ff versus Aa is shown for three different values of the ratio 
Q n /uo. As shown in (17) and (18), while the classical limit is the same for all the situations, 
different behaviours appear in the quantum regime. The minimum effective uncertainty is 
maximum at the resonance condition Q n = uj. In this case it is always Aa e // > ^jfi/muj. 
When Q n ^ u the minimum effective uncertainty estimated by the intersection between 
quantum and classical limits decreases as (f2 2 — uo 2 )^ 1 . 

The comparison between analytical and numerical results for Aa e // are shown in Fig. 3 
for two different choices of the measurement parameters. Also shown are classical and 
quantum behaviours, corresponding to the limits expressed in (17) and (18). The difference 
between numerical and analytical results is less than 0.1% and it can be further reduced by 
choosing higher resolution space-time lattices. 

The numerical accuracy estimated above allows us to perform meaningful computations 
for a nonlinear oscillator represented by the Lagrangian 

L = —x 2 — ^* 2 - U (19) 
2 2 4 K J 

For comparison with the previous linear case we have chosen the initial and final states of 
the form (13) and a measurement output of the form (14). In this case we expect both non- 
Gaussian behaviours for the wave function at t — r and the propagator. This also means 
that the distribution P(e) is not a Gaussian function. An equivalent width for P(e) may be 
introduced through the definition 



Aa <» s 7mo) J- mdt (20) 

The computed Aa e // versus Aa are shown in Fig. 4 for two different values of the nonlin- 
earity coefficient (3. The comparison with the linear situation having the same parameters is 
also shown. It appears that the effect of the nonlinear term is in the direction to enlarge the 
region in which the classical approximation is meaningful. The quartic term concentrates 
the final wave function near the result of the measurement and this implies that the effect 
of quantum noise, toward a spreading of the most probable paths, is reduced. 

In all the previous considerations we have chosen to deal with Aa which is time inde- 
pendent. A more general class of continuous measurements is obtained by considering time 
dependence for Aa. In particular to recover chaotic dynamics through quantum measure- 
ments a particular conditions on the kind of measurements must be satisfied, namely the 
measurement process has to be a quantum nondemolition process for the observable under 
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monitoring [4]. Quantum nondemolition strategies for linear oscillators have been already 
analyzed in the framework of the path integral with continuous measurements [8]. The ap- 
plication of the technique described here for quantum nondemolition measurement processes 
on nonlinear systems will be the subject of future investigations. 

Acknowledgements: One of us (M.B.M.) is indebted to Prof. G. Immirzi and Prof. F. 
Marchesoni for the kind hospitality in the University of Perugia where a part of this work 
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FIGURES 
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FIG. 1. A continuous measurement of position leads to a corridor centered at the measurement 
result a(t) where the most probable paths lie. 
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FIG. 2. Effective uncertainty versus instrumental uncertainty for three different values of 
the ratio Q, n /oj in the case of a linear oscillator. The dashed line is the classical limit Aa c , the 
dot-dashed lines are the quantum limits Aa q . The value of the other parameters is also indicated. 
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FIG. 3. Numerical (dots) and analytical (solid line) effective uncertainty versus instrumen- 
tal uncertainty for a linear oscillator. Asymptotic classical (dashed) and quantum (dot-dashed) 
behaviours are also shown. 
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FIG. 4. Numerical (dots) effective uncertainty versus instrumental uncertainty for the non- 
linear oscillator at two different values of (5. The solid line is the behaviour of the corresponding 
linear oscillator obtained for (3 = 0. The dashed line represents the classical limit. 
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